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We extend field theoretic variational perturbation theory by self-similar approximation theory, 
which greatly accelerates convergence. This is illustrated by re-calculating the critical exponents of 
0(Af)-symmetric (f)'^ theory. From only three-loop perturbation expansions in 4 — e dimensions we 
obtain analytic results for the exponents, with practically the same accuracy as those derived recently 
from ordinary field-theoretic variational perturbational theory to seventh order. In particular, the 
theory explains the best-measured exponent a ft; —0.0127 of the specific heat peak in superfiuid 
helium, found in a satellite experiment with a temperature resolution of nanoKelvin. In addition, 
our analytic expressions reproduce also the exactly known large-A'^ behaviour of the exponents u 
and 7 = !/(2 — 77) with high precision. 

PACS numbers: 74.20.De,74.25.Fy,74.40.-|-k 

I. INTRODUCTION 

The precise calculation of critical exponents of phase transitions is an important theoretical task. On the one hand, 
these exponents provide us with basic information on the behaviour of thermodynamic quantities in the vicinity 
of critical points. On the other hand, such calculations require the development of new mathematical techniques 
to master the resummation problem of divergent perturbation expansions (see e.g. The comparison of the 

calculated exponents with experiment serves as a test for the validity and accuracy of the mathematical methods. 

Recently, apowerful method has been developed by one of the authors called field-theoretic variational perturbation 
theory which converts divergent weak-coupling into convergent strong-coupling expansions. This method 

presents a more powerful alternative to the previously used method of Pade-Borel resummation. The higher power 
stems from an explicit knowledge of the convergence behavior expressed in terms of the Wegner exponents lu of 
approach to the strong-coupling limit. The higher accuracy has been amply demonstrated by calculating the critical 
exponents in particular, by predicting the most-accurately known exponent a « —0.0127 of the specific 

heat peak in superfiuid helium, found in a satellite experiment with a temperature resolution of nanoKelvin. 

An important feature of any resummation method is the convergence of the renormalized sequence of approximants. 
This can be studied by considering the analytic properties of the sought function with respect to coupling and by 
involving the corresponding dispersion relations [ij, |3 • The field-theoretic variational perturbation theory was shown 
to possess an exponentially fast convergence (see the detailed proof in Q). 

In an independent development, the other of the authors has developed a completely different general resummation 
scheme called self-similar approximation theory [^ linLlllIll2tlT^ . Also this method exhibits a fast convergence, which 
has been demonstrated for a variety of problems in quantum mechanics, statistical physics, and mathematical finance 
(see review- type papers 0,^3)- 

The aim of the present joint paper is to combine the two approaches. The combination is expected to have the 
fastest convergence so far. Indeed, we shall find, that the combined method is capable of reproducing from three- 
loop expansions the above mentioned seven-loop result of field-theoretic variational perturbation theory for critical 
exponents, including the most sensitive critical exponent a. 

In Section II, we give a brief reminder of the basic formulas of field-theoretic variational perturbation theory, which 
will be used as a basis for a further acceleration of the convergence via self-similar approximation theory to be reviewed 
in Section III. In Section IV we develop the combination of the two methods, which is then applied in Section V to 
calculate the critical exponents of the 0(A^)-symmetric field theory. 

II. FROM WEAK TO STRONG COUPLING 

Physical quantities of interest are usually derived from theories as divergent series in powers of some bare coupling 
constant gs- These provide us with reliable results only for very small gs- Critical phenomena, however, take place at 
infinitely large in comparison with the mass, the inverse length scale of the fiuctuations [ij . In order to overcome 
this difficulty one has to reorganize the divergent weak-coupling series into a convergent strong-coupling expansion. 
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Such a reorganization is provided by the field-theoretic variational perturbation theory 0, 0, 0, 13 1 briefly summarized 
in this section to recall the principal formulas needed in what follows. 

Consider a real function /{gs) of a real gB, whose limit /(oo) we want to find from a divergent weak-coupling 
expansion up to order L, 

L 

f^''\9B) = Y.f^^9l (93^0), (1) 

n=0 

with L = 1,2,3,... enumerating the maximally available order. Our aim is to find the behaviour of figs) at gs oo. 
The expansion coefficients, grow factorially with n so that the series could have sense only for very small gs- 

Field-theoretic variational perturbation theory is based on the introduction in of a variational parameter K by 
the identical replacement 

9B „ l-K^ 

[K-^ + gBr)i gs 

where g is a parameter related to the critical Wegner exponents cj, which in the rcnormalization group approach to 
critical phenomena governs approach to scaling. The parameter q in this paper corresponds to q/2 in the original 
work Isj. After the replacement, the series 1^ is re-expanded in powers of gs at fixed r, and at the end r is again 
replaced by (1 — K^)/gB- This procedure introduces an artificial dependence on the dummy parameter K which is 
fixed by searching for a plateau in K which becomes fiatter and flatter for increasing order. The plateau is horizontal 
only for the correct choice of q, and this condition will determine the Wegner coefficient lo 

In the upcoming calculations, we shall work with a slightly different but completely equivalent replacement 

1 - gBr)i gB 



where a is defined as the function oi gB'- 



a ^ a{gB) ^ 1 - ( — ) . (4) 



s 

9B 



Following the rules of field theoretic variational perturbation theory, we have to form the functions 



- gBrY 



(5) 



to be calculated with the prescription that the terms g'% = s"/ (1 — gBi')^'' in the truncated series (Q) are re-expanded 
systematically in powers oi gB up to 5^ After this, we replace again r <y I gB and optimize the resulting function 
in the variational parameter s. Using the binomial expansion 



L — n 



r(m + i)r(p- m + 1) 

we obtain explicitly 

F^^\gB.s.,q) = J2J2 (-^)'" ■ (7) 



L L- 



n— m—O 



This must be optimized in s, yielding order-dependent function s^^\gB) and an associated cr^^^(gs) = 1 ~ 

{s^'^\9B)lgBf'. 

Note that Eqs. (2) and (3) are identities and not directly related to the functions appearing in the Symanzik-type 
transformations |l7l | , in spite of a certain similarity. 

Our aim is to find the behaviour of Eq. (U)) in the strong-coupling limit gs —> 00. From Eqs. (PJ we observe 
that a ^ I as gs ^ 00 since, as we shall see, the optimal s is finite. This allows us to calculate from Eq. Q as 
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finite approximant F*-^^ (oo, s, g). In what follows, we limit ourselves to the approximants of third order, since then 
all calculations can be done analytically. The first three approximants are explicitly 

F(i)(cx3,s,g) = /0 + /1 s, (8) 
F(2)(oo,5,r7) = F<^'\^,s,q)+q hs + f2S^ , (9) 

F(3)(oo,s,<7) = F'^^'>{^,s,q) + ^q{l + q) hs + 2qhs^ + hs^ . (10) 

The optimal values s*-^' = s'^^^(oo) are found either by extremization 

d 



= 0, (11) 



S — S 



or, when the latter has no real solutions, from the turning points 

= 0. (12) 
To second order, there exists an extremum at 

s^'^=-{l + q)^- (13) 

To third order, there can be two possibilities. There is an optimal extremum s — s*-'^^ at one of the roots of the cubic 
equation 

3/3 + 2(1 + 2(7) /2 s + i (1 + q)i2 + q) f, ^ , (14) 



and an turning point at 



5(3) = -(l + 2g) A. (15) 



Usually, conditions 111|) and (|12|l yield optimal values of s^^'^ alternatively for odd and even orders L, respectively 
(see 1311 mi Si, 18]), and this will be the case in the upcoming applications of this paper. 
After determining s^^\ we obtain the optimized approximants 

F^^^°P\oo,q) = F^^\oo,s'^^\q) . (16) 

To second order this is 

F(2)°Pt(oo,g) = /o-(l + -?)2-^, (17) 

4/2 



and to third order, with s^^^ of Eq. (QSJ (see Ref. i^). 



opt(oo, g) = /o - (1 + q){l + 2g)(2 + g) M + (1 + 2q)3 ^ . (18) 

0/3 ^ ' J3 



For each approximant H16() . we must also specify the parameter q. If the Wegner exponent were known from other 
sources, we could use this. Otherwise we must determine it order by order, which yields an L-dependent result 
q = q''^\ so that the final approximants will be 

y(L)opt ^ ^(L)opt^^^^(L)) _ (19) 

The determination of g^^^^ proceeds as follows. If we expect the function figs) to be finite in the strong-coupling 
limit, which is the case for the critical exponents, then the logarithmic derivative, to be referred to as /3-function Q], 

P{9b) ^ '-^If^ (20) 
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must tend to zero for gs ^ oo: /?(oo) = 0. From the expansion it is straightforward to derive 

P^''\9B) = Y.P-9l- (21) 

Tl = 

Depending on whether /g is nonzero or zero, the coefficients /3„, up to third order, are given either by 



Jo Jo Jo 



f3o = 0, /3i = ^ , (32 = 2^ - ^, (22) 



/?3 = 7T-3^+3f (/oy^O), (23) 
Jo Jo JO 



or by equations 
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/3o = 1 , /3i = ^ , /32 = 2 ^ - ^ , (24) 
/i Ji Ji 

P3 = 4 -3^+3^ (/o = 0) . (25) 



/i^ ft /: 



Now we treat the expansions P^-'^^gs) in the same way as before /^^Hss)- We form the approximants S*^^^ °p*(oo, q) 
similar to the way of deriving Eq. H16() . This is set i?'^^-' °p*(oo, 5) equal to zero to ensure /3(oo) = in each 
approximation. This determines the proper parameters q = g'^^. For instance, 



1^'' - 2 - 1 . (26) 

Note that the logarithmic derivative for determining q^^^ can be formed from any function of gs with a constant 
strong-coupling limit l^J , i.e., from any critical exponent, not just from the function figs) we want to resum at the 
moment. Usually, the function gR^gs) relating gB to the renormalized coupling constant ga is most convenient, since 
it is known to highest order. 



III. SELF-SIMILAR APPROXIMATION THEORY 

Self-similar approximation theory [^ ll0lllllll2llT^ is based on constructing a sequence of optimized approximants, 
which contain instead of a variational parameter a trial function. The general idea of deri ving convergent sequences 
of optimized approximants with the help of trial control functions has been suggested in |2fl|. The first step in the 
optimization procedure is reminiscent of the Euler-Lagrange variational method. But while the latter is a single-step 
procedure [2l|, the optimized perturbation theory runs via a sequence of better and better approximants. 

In the last section we have shown how to calculate a sequence of trial functions {-F^^'((7_b, s, q)} by field-theoretic 
variational perturbation theory [see Eq. iTt Ii]. F rom these s and q can be determined as functions oigs by optimization. 
In self-similar approximation theory |9|, llOl IllL Il2l |l3j , the approximants of different order are considered as a flow on 
the manifold of approximants, in which order L of the approximation plays the role of a discretized pseudotime. In 
this interpretation, the sequence of approximations behaves like a dynamical system. The higher approximations will 
be obtained by improving the entire optimized functions s^^'' (55)1 even if we are only interested in the strong-coupling 
value /(oo), for which the previous method required only an optimal parameter s''^)(oo). Thus we have to perform 
the optimization procedure for all gB before going to the limit gs 00. Instead of Eqs. and H12f) . we have to 
solve the full extremality condition 



= 0. (27) 



and, if this has no real solution, the turning point condition 

92 



F^''H9B,S,q) 



= 0, (28) 
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to find the lowest approximation for the trial functions s^^-' = s'^^^gs)- More explicitly, we could also record the 
parameter q at which the optimization is done in the arguments and write the solution as s^^\gB,q)- But we shall 
refrain from doing so to avoid cluttering the notation. For the same reason we shall omit, for a while, the argument 
q in F(^)(oo, s, q). 

Starting from the trial functions s^^'(gB), we construct an approximation following the general scheme developed 
in H [IS m [Hill El El. We define the reonomic function gs = g^^\4>) by the reonomic constraint 

J^^'HffB, 5(^^35)) -0, (29) 
where F^^^ is the lowest nontrivial function in the sequence {i^*-^^}. Now further define an entire sequence of functions 

y(^)(0) ^ F'^^\g'^^\<j^),s^^\g'^^\^)) , (30) 

with the initial term y^^^cj)) = (j). The set of all functions y^^\4>) for L = 1, 2, 3, . . . constitutes a space y C y'^^^cj), q) 
called approximation space. The pseudotime evolution in this space forms a group of self-similarity transformations: 

/^+P)(0) = y(^)(/^')(0)) . (31) 

The property of self-similarity (|31|l guarantees the existence of a fixed point y* = y^^-** [l^ll5| . which has the property 

r = y^^^y*) . (32) 

More explicitly, the fixed point satisfies 

y(L). ^ F^Lh(^giL)^ ^ s(^)(g(,^')) , L > 1. (33) 
It defines the desired self-similar approximant 

/(^)*(gB) ^ F(^)*(5b(0),s(^)(5b(0))) . (34) 
In order to find the fixed point we define a pseudovelocity of the approximation sequence by the finite difference 

t;(^)(</.) ^F^+i(.gl,^) (</.), s(^)(0))-F(^) (3(3^) ((/.), sW(0)) . (35) 

If the {y'^^^ with discrete L = 0, 1, 2, . . . were a flow function {y^*-*} of a continuous time t > 0, it would follow a time 
evolution equation 

|2/W(0)=«(*)(y(*)(0)) . (36) 



The integral form of the latter can be presented as the evolution integral 

(37) 



ly(L) v(^){(f>,q) L ' 

If the parameter q is unknown, it must be determined from a simultaneous treatment of the /3-function (|20|) . In this 
case we determine a sequence of optimal parameters q — g*-^^*, leading to the self-similar approximants 

/(^)*(.gs) ^ F(^)*(gB, s(^)(gB), 9^"^^*) • (38) 

For the purpose of determining critical exponents we are only interested in /(gs) at gs — > 00 and go to the limit of 
Eq. ||2HI) yielding 

/(^)*= hm /(^)*(.9b). (39) 

The self-similar approximant (|39|l replaces the previous optimized approximant p9l) of field-theoretic variational 
perturbation theory. 

From the definition of the pesudovelocity (|35|l it follows that for the calculation of the i-order self-similar approx- 
imant Z^^-**, we need to know L + 1 orders of the expansion in (^). If L is the last available order, we shall use as an 
L + 1st approximation the average of the previous ones: 

j(L+l)* _ 1 j-j(L-l)* ^ . (40) 

This approximation is expected to be reliable if the approximants tend to the limit L ^ cx) in an alternating fashion, 
once from above and once from below. This is not a priori ensured, but happens in many examples (but there exist 
also counter-examples). In the series for the critical exponents to be treated here this seems to be true. Thus Eq. 
(|40|l will be used to get the highest approximant. More refined mathematical foundation for the usage of Eq. 139|) is 
given in Refs. [2^l23|. 



6 



IV. COMBINING SELF-SIMILAR AND VARIATIONAL THEORIES 

Let us now be explicit and improve the convergence of the sequence {J(^)°p*} of variational perturbation theory 
derived in H16|l - (|19|1 by self-similar approximation theory to obtain a new sequence {Z^^-**}. The improvement is most 
drastic at the initial stages of the procedure, when i < 3, so that we shall restrict ourselves to these low orders. An 
additional advantage is that all formulas up to third order can be derived analytically. 

Recall that, in contrast to Section II, we do not consider from the beginning the limit of gs oo, but retain the 
full gs-dependence of the functions Q : 

F^^\gB,s,q) = /0 + /1S, (41) 

= F'^'\gB,s,q)+qhas + f2s'', (42) 

F^^\gB,s,q) = F^^\gB,s,q) + ^q{l + q)ha\s + 2qf2as^ + hs\ (43) 

which reduce to Eqs. |(S||- H1U|) for gs ~> 00 since then cr — > 1. For arbitrary gs we must optimize i^(^) in s. Since a 
depends on s via the relation I^J , we may look for the extremum in the two s and a while satisfying the condition 

da a — I , ^ 

— = . 44 

as qs 

If this is done with the function F^^\ we obtain an optimal function s'^^'^(gB)- From this we calculate the approximant 

F'^'^^°^\gB,q) ^ F^^\9B,s^''\gB),q) • (45) 

To lowest order L ^ 1, an optimal function dos usually not exist. In this case we shall use the next higher existing 
s^^Hfs) to define the lowest approximant. In principle we could, of course, form an entire off-diagonal matrix of 
variational functions 

F(^^'^'\gB,q)^F^^\gB,s^'^'H9B),q) • (46) 
of which the functions l|16|l in Section II arc only diagonal elements 

F(^)°p'(oo, q) = F(^'^)(c50, q) . (47) 
The optimal function 5(^^(5^) is determined by the extremality condition H27|) . which amounts to the equation 

(l + 9)/ifT(5s) + 2/2S = 0. (48) 
From this we obtain the variational expression 

F(^^°P*(<?s,g) = F^'-^\gB,q) = fo + h .s(^) + ^ /2 ^^^'^ ■ (49) 

In addition, we determine a lowest approximant with a first-order trial function s'-^-'((7s) = s^'^Hqb), 

F^''^°'P'{9B,q)^F(''^HgB,q) = h + fiS^'^ . (50) 
The third-order trial function s^^\gB) is given either by the extremum-point condition (|27|l . which yields 

i (1 + q){2 + q)ha^ + 2(1 + 2<z)/2a,s(3) + 3/3s(3)2 = , (51) 

or, when Eq. (|51|l has no real solution, by the turning point condition (|28|l . by solving 

(1 + g)(2 + q)ha{a - I) + 2(1 + 2g)/2 s^^'^a - 1 + ga) + 6/395^'^' = . (52) 

We now construct the approximation cascade following the rules of Section III. Accordingly, we set y'^^^ = F'^^\ 
y{2) — ^(2)^ and — F^^) for L > 2. Therefore the reonomic-constraint H29|l reads fo + fis = 0, which defines the 
reonomic function ''('/') through the relation 

s^'\9Bm = ^ . (53) 
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The first-order psudovelocity as defined in Eq. H35|l is u^^^ = i^(2,i) _ f'^^'^\ which reads explicitly 

i;W(0,g) = Ai(0-/o)2, (54) 

with 

(1 - g)/2 

Using the evolution integral (|37(l . we find from self-similar approximant H33(l the expression 



^1 ^ • (55) 



where F^^^^) 

is given by Eq. H49() . In the strong-coupling limit — > cx), this reduces to 

(5 - qyf2 



(57) 



Let us now determine the value of the parameter q, following the procedure described at the end of Section II, but 
with the difference that now we construct the self-similar approximants for the expansion (|21|l of the /3-function, which 
will be denoted by B'^^^*{gB, q)- These have the same form as F^^^*{gB, q), except that the expansion coefhcients /„ 
are replaced by /3„. The parameters g^^-** are defined by the boundary condition 

B^^>{oo,q'-^>) =0 . (58) 

To first order, the result is 

Substituting this into Eq. we obtain f^^^*{gB) = F^^''* {gs, q^^'^*), in agreement with Ea. (|Mll . The final result is 
given by Eq. that is, by the value /(i)* = /(i)*(oo). 

The second-order velocity (p?5| is w^^^ — i^(3,2) _ p{2,2) ^ where 

F^'''\gs,q) = F^'''Hgs,q) + (^ + f ff^ -^lf^ ,(2)3 (6o) 

U + q)ji 

which is valid for any gB- In particular, in the limit gs oo, where 

A3,2)f^^ „\ _ f /T , „^3 fi f /1/3 



The explicit form of the velocity u^^^ is now 

v^-^\<l>,q)^A,{<l>~fo)\ (62) 

where 

, _ (1 + q)fih - 2qfi 

= il + q)ft ■ ^''^ 

From the evolution integral H37|l we find here 

F'^'^igs, q)^fo+ F(^^'\gB,q)-fo 

^l-A2[Fi^^^K9B,q)-fo? 

The value g'^^^* follows from the strong-coupling condition H58|) . which leads to the equation: 

(1 + q^'>)\0oPlP3 - WlPl + Pt) + 2(1 + q^'^^* fPlPl - KSplPl = . (65) 

Inserting the appropriate solution g'-^^* into Eq. (|^ . we obtain f'-^^^*{gB) = F^^^*{gB,q''^''*)- And the limiting value 
is = /(2)*(oo). If only three orders of the expansion (QJ are available, then /f^)* jg defined by Eq. H4l)|) . as 
explained at the end of Section III. 
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V. APPLICATION TO CRITICAL EXPONENTS 

The above theory will now be applied to evaluate the divergent perturbation expansions of the critical exponents of 
0(A^)-symmetric 0'*-field in 4 — e dimensions. The expansions are power series in the bare coupling parameter gs/u^, 
where /i is some mass parameter to make gs/l^'^ dimensionless. They can be found up to six loops in the textbook [jj]. 
In this field-theoretic context, the parameter q in the transformation (j^J is directly related to the Wegner exponent 
uj |24l | which characterizes the strong-coupling behaviour of the renormalized coupling 

9(93) ^ g(oo) - const X -— {93/ fJ^' 00) . (66) 

9b 

The relation is 

u^e/q. (67) 

In the sequel, we shall set fi — 1. Starting point is the expansion for the renormalized coupling constant which we 
shall limit to g%, for simplicity; 

9R{gB)^gB + C2g% + C3g% + C4g% , (68) 
where gs ^ and the coefficients are 

N + 8 + 3iV+14 

Co = , Ci = 1 , C2 = — , C3 = — \ , 

3e 9e^ be 

(iV + 8)3 4(A^ + 8)(3iV+ 14) SSA^^ ^ 922Ar + 2960 + 24(27iV + 88)C(3) 
^^^"^7^5 ^2 ' (69) 

with C(z) being the Riemann zeta function. The logarithmic derivative of (|68|l yields the /3-function 

^(3) (^^) ^ dl j9i9B) ^ ^ ^ ^^^^ ^ p^g, ^ ^^g, ^ (70) 
d log 55 

with the coefficients 

/?! = C2 , 132 = 2C3 - , = C2 - 3C2C3 -|- 3c4 . (71) 

In second-order variational perturbation theory, we have q^^^ given by Eq. H26f) . which yields 

=2yTT^-l, (72) 

where the notation 



3(37V-|- 14) 
{N + 8)2 



P = , o^2 (^3) 



is used. The first-order self-similar approximant q'^^'^* is defined in Eq. H59|l . resulting in 



^(D* ^ V(l+pe)(9 + 5p£)-l _ 



The corresponding Wegner exponents are 



,(2) = 1 



2VTTpi- 1 

in the variational perturbation theory, and 



(75) 



V(l+pe)(9-t-5pe)-l 
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in the self-similar approximation theory. 
The second-order optimized approximant 



^ ^ 4(7V + 8) ' ^ ' 



corresponding to the renormalized coupling (|68|l . with the parameter (|72|l . becomes 



^(.)opt^3(^^ (78) 



N + 8 

according to Eq. (|19|1 . The first-order self-similar approximant 

with q^^'>* from Eq. ifTH) . results in gl ~ Gl{oo,q^^^*), as in Eq. It turns out that, because of the equality 

(l+g(l)*)2 



5- 



1+pe, (80) 



the values g^^^* and g^^)"?* coincide. However, f/^^^*, following from Eq. is different from (;(3)opt_ 

We now turn to the perturbation expansions of the critical exponents v and 7. Other exponents need not be treated 
since they can be found from the above using well-known scaling relations (see e.g. 0,11^11^). We begin with v^^ 
for which we use the expansion up to 



V 



-1(3) 



with the coefficients 



/o + /i.9B + hg\ + hg\ , (81) 



h = - 



iV + 2 
108 



4(7V + 8)2 2(19iV+122) 



3(5Ar + 37) 



Following the above procedure, we get the optimized strong-coupling value 

(A^-H2)e 



^-1(2)* ^2-(l + g)2 



2(2iV -I- 16 5e) 



(83) 



(84) 



which, for cp''^ from Eq. H72|l . gives the variational perturbation result v 1(2) op* — y °p*(oo, g*^^)). Its self-similar 
improvement reads 



2(l + (7)2(Ar + 2)6 



(5-g2)(2iv+16 + 5e) 
After inserting q^^^* from Eq. 1)74(1 . we obtain 

2(iV + 2)(l+pe)e 
2(iV + 8) + 5e 



(85) 



(86) 



Again, it turns out that v'^^^'^* = jy-i(2)opt^ ^mt i/-^^^)* ^ j,(3)opt^ .^^j^^j^ i/-i(2)* defined by Eq. 

Finally, we resum the perturbation expansion for the critical exponent 7 = ^{2 — rf), which reads in the form P|: 



with the coefficients 



/o = 1 , /i = 



7(35) =i /o + fWB + ./2.g| + /3.9|j , 

2(iV-H8) 



(87) 



N + 2 
6 



/2 = - 



36 



4-iV 
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73 - 



8{N + 8y 4(106 + iV - 2Ar2) 



432 

Here we find tlie optimized approximant 



194 + 7V(2Ar + 17) 



(89) 



and the self-similar approximant 



4[2(7V + 8) + (4 - iV)e] ^^^^ 



^W* ^ 1 + + (91) 
^ 2{N + 8) + - N)e ■ ^ ' 



Again 7(1)* coincides with the variational perturbation result 7(2) opt _ ^(2)optf^^^ g'-^-'), whereas 7^^)* does not equal 

^(3)opt_ 

By construction, the self-similar approximants /(^)* obtained from the evolution integral H30|) possess the same 
e-expansion, up to the given order L, as the optimized approximant J(^)°p* of variational perturbation theory. This is 
evident^rom expressions (|56() and H64() . For the variational perturbation results, on the other hand, it was shown in 
Ref. jj, y, llSl that all expansions in powers of e coincide with the expansions derived in the renormalization-group 
approach to critical phenomena |l|. As a consequence, also the presently derived self-similar approximants /(^)* 
possess the exact e-expansions. This can easily be verified by an explicit calculation. 

It is interesting to compare the l/A^-expansions with the self-similar approximants for larger N. These expansions 
for the critical exponents tu, v, 7, and ry are presented in Fig. 1, where they are compared with our self-similar 
approximants as well as with the results of the sixth-order variational perturbation theory |l|, and with those of 
Pade-Borel resummations j^^,!!^. Our third-order results have the same accuracy as those of sixth or higher orders, 
obtained by other resummation techniques. In the limit N 00, our exponents coincide with the known exact values 

D-A 1 2 D + 2 , , 



„=-^—^, ,y = 0, u; = 4-D, (93) 

where D is dimensionality. 

Critical exponents for finite N have been calculated by Pade-Borel resummation methods based on six- and seven- 
loop expansions in D = 3 dimensions 29, 30, 31. 3^ or in five-loop expansions in = 4 — e dimensions 33, 34, 35, 3^. 
Different variants of the optimized perturbation theory have been used HE HE QE| ■ Self-similar exponential 
approximants were given in |4]| . Computer simulations, based on the Monte Carlo lattice studies were presented in 
[421 . The available results have been reviewed in Rcfs. 1, 43. 44]. 

A list of our results from the third-order self-similar improvements of variational perturbation theory for the critical 
exponents is given in Table 1. The exponents v, 7, and lo are calculated directly from their series, as is explained in 
the text. The other listed exponents are obtained from the scaling relations 

a = 2-iyD, /3 = ^ (Z? - 2 + 77) , ^ = ^(2-7^). (94) 

The error bars are defined by the difference between /(■^)* and /(-^^i)*. Our results are compared with those of 
the field-theoretic variational perturbation theory based on six-loop (for > 3) and seven-loop (for A^ = 0, 1, 2, 3) 
expansions calculated in three dimensions in Refs. 0, S IBlJ and listed in the book We also show the exponents 



recently obtained by Pade-Borel resu mmation l28l. as well as earlier results |29H3nH3lLl32ll33ll34l35ll3fi| . base d on six - 

Dns Liii [ 



loop expansions in D ~ 3 dimensions '2^,'3^,'3l','32'| and on five- loop expansions in D = 4— e dimensions j33Ll33. l35'.'3^ 

The numbers in parentheses indicate the highest calculated approximation (seventh order for n = 0, 1, 2, 3 and sixth 
order for A^ > 3), from which the effective extrapolations to infinite order were obtained as described in the book jj]. 
Comparing the results, we see that our third-order self-similar approximants yield practically the same values for the 
critical exponents as the values derived by other resummation techniques of sixth or seventh order. In the limiting 
cases of A^ = — 2 and N — 00, our results coincide with the known exact values of the critical exponents. 

The critical exponents obtained by our new method from three-loop perturbation expansions are all in good agree- 
ment with all experiments. Unfortunately, most of them are not sufficiently accurate to distinguish between different 
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theoretical approaches. The most accurately known experiment is the measurement of the specific heat of hquid hehum 
with nanoKelvin temperature resolution near the lambda point, which were performed in a satellite orbiting around 
the Earth 45, 4(|. The specific-heat exponent initially extracted from the data in was a — — 0.01056 ± 0.0004. 
This differed slightly from the result of seven-loop variational perturbation theory obtained from three-dimensional 
(/)^-theory which yielded a — —0.0129 ± 0.0006 0. However, a recently performed re-analysis of the data in Ref. [i^ 
found a — — 0.0127 ± 0.0003, thus confirming with great precision the theoretical result of Ref. 1^. Our present result 
a ~ —0.0124 obtained in third-order self-similar-improved variational perturbation theory is again in perfect agree- 
ment with the latest experimental result. This is quite remarkable since the five-loop calculations in 4 — e dimensions 
gave the value a — —0.013 This illustrates the acceleration of the convergence of variational perturbation theory 
by the self-similar improvement developed in this paper. 

In conclusion, we have developed a novel method for resuming divergent perturbation expansions. It combines 
field-theoretical variational perturbation theory with the self-similar approximation theory, and accelerates greatly 
the convergence of either method alone. The acceleration is especially useful if only low-order expansions are available 
due to the complexity of the problem. Up to third order, all results are found analytically. The method was illustrated 
by calculating the critical exponents whose third-order approximants are found to be practically the same as the sixth- 
or seventh-order approximants of other resummation techniques. The specific-heat critical exponent a for = 2 is 
found to be in perfect agreement with the most accurately measured experimental value of a for superfiuid helium. 
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FIG. 1: Solid curves show our fourth-order approximations to u, i/, 7, rj. Short-dashed is second-, long-dashed curve is third- 
order approximation. Thin dots show sixth-order approximation of the textbook fat dots the extrapolations to infinite 
order. The dash-dotted lines in the second and third figures are interpolations to the Pade-Borel resummations of 
(where ui was not calculated). Their data for r/ scatter too much to be represented in this way-they are indicated by small 
circles in the fourth figure. The dotted curves show 1/A'^ -expansions of all four quantities. Note that our results lie closer to 
these than those of S.A. Antonenko and A.I. Sokolov. The solid r;-curve was calculated in the textbook Q (see Fig. 20.2). The 
exact large- limits are ujn=oo = 4 — D, vn=oo ~ 1/{D — 2), 7]v=oo — 2/{D ~ 2), and i]m=oc = 0. The exact values at A = —2 
are vn=-2 = 1/2, 7]v=-2 = 1, and riM^-2 ~ for all D. 



